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Abstract. We explore the quantum dynamics of heteronuclear atomic collisions in 
waveguides and demonstrate the existence of a novel mechanism for the resonant 
formation of polar molecules. The molecular formation probabilities can be tuned 
by changing the trap frequencies characterizing the transverse modes of the atomic 
species. The origin of this effect is the confinement-induced mixing of the relative 
and center of mass motions in the atomic collision process leading to a coupling of 
the diatomic continuum to center of mass excited molecular states in closed transverse 
channels. 
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1. Introduction 

Dimensionality plays a crucial role for interacting ultracold quantum gases [H El [3l H] . 
For a strong two-dimensional confinement, i.e. a tight waveguide, and a close to resonant 
scattering process the atoms can feel the trapping potential within an individual atom- 
atom scattering process. A striking manifestation herefore is the confinement-induced 
resonance (CIR) effect [5j. The narrow transversal confinement leads to a singularity of 
the coupling constant for the corresponding effective longitudinal motion and as a result 
we encounter a ID strongly interacting gas of impenetrable bosons that can be mapped 
at resonance on a gas of free fermions, the so-called Tonks- Girardeau gas [6]. The CIR 
has been interpreted as a Feshbach resonance involving the atom-atom continuum and a 
bound state of a transversally excited channel of the harmonic waveguide [3 [HI [9]. CIRs 
have subsequently been discovered for three-body [101 E] ^^^^ four-body [12] scattering 
under confinement and for a pure p-wave scattering of fermions [13] . The CIR behaviour 
has been found experimentally for s-wave scattering bosons [TH [15] and for p-wave 
interacting fermions [16]. Recently a dual-CIR, which leads to a confinement- induced 
transparency effect [171 [H|, has been demonstrated. It is characterized by a complete 
suppression of strong s- and p-wave heteronuclear atomic scattering in 3D because of 
the presence of the waveguide. 

It is to be expected that waveguides with an anharmonic and/or anisotropic 
transversal confinement will lead to novel ultracold atomic collision properties and 
consequently to an intriguing many-body dynamics of the quantum gas. Very first 
results in this direction include the observation of a modified CIR due to a coupling of 
the center of mass (CM) and relative motion of the atoms [T8ll8]. This nonseparability 
can be either achieved by employing anharmonic waveguides or, in case of heteronuclear 
collisions, by the fact that the two species experience different harmonic frequencies. 

In the present work we show that the quantum dynamics of the coupled CM 
and relative motion for two-atomic species feeling different confinement potentials, i.e. 
waveguides, exhibits a resonant formation process of ultracold molecules. The latter 
is analyzed in detail and its existence is found to be independent of the specifics of 
the atomic interaction. As a consequence it is anticipated that heteronuclear atomic- 
molecular ID quantum gases in an anisotropic waveguide should explore a rich many- 
body quantum dynamics. 

2. Wave-packet propagation method 

We employ a wave-packet dynamical approach to study the atomic scattering in 
waveguides which has been developed recently [ITJ [8]. The collisional dynamics of 
two (distinguishable) atoms with coordinates ri,r2 and masses mi,m2 moving in the 
harmonic waveguide with the transverse potential | X^i ^i^fPi {pi = TiSinOi) is described 



Quantum dynamics of resonant molecule formation in waveguides 



3 



by the 4D time-dependent Schrodinger equation {h = 1) 

d 

«^^(P/?, r, t) = H{pR, r)ij{pn, r, t) 

with the transformed Hamiltonian 

H{pR, r) = H^{pr) + H^{r) + W{pn, r) + V{r) . (1) 

Here 
and 

Ipor^ Ipv^ 2 mi m2 
describe the CM and relative atomic motions, V{r) describes the atom-atom interaction, 
Pr and r = Ti — r2 (r, ^, 0) i-^ (p, 0, 2;) are the polar radial CM and the relative 
coordinates and M = mi + m2, p = mim2/M. The transformed Hamiltonian ([1]) is 
given in a rotated frame which exploits the conservation of the component (Li + L2) of 
the total angular momentum and allows to eliminate the angular dependence of (f)R of 
the CM degrees of freedom P [19]. Our investigation will focus on the case M^^ = 0, 
the later being the quantum number belonging to j-^^- 
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The term 
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W{pR, r) = p{uJi - uj^ypR sin 9 cos (4) 

in the Hamiltonian ([1]) leads for two distinguishable atoms that feel different confining 
frequencies uji ^ uj2 io a, nonseparability of the CM and relative atomic motion. We 
integrate the Schrodinger equation from time t = to the asymptotic region t +00 
with the initial wave-packet 



ip{pR, r,t = 0) = Nr^/p^exp{-^^ ~ 2a^ 2a? — ^ 

representing two different noninteracting atoms in the transversal ground state of the 
waveguide with = (l/mjtUj)^/^ and the overall normalization constant N defined by 
< il:{0)\ip{0) >= 1. We choose Zq —00 to be far from the origin z = and — 00 to 
obtain a narrow width in momentum and energy space for the initial wave-packet. Our 
wave-packet moves with a positive interatomic velocity vq = k^/ p = ^2e\\l p thereby 
approaching the scattering region located at 2; = 0. In the course of the collision, 
the wave-packet splits up into two parts moving in opposite directions z ±cx3. 
We model the interatomic interaction Vir) via the Lennard- Jones 6-12 potential 
Vir) = Cu/r^"^ — Cq/t^. Note, that the effects discussed in the following exist for even 
qualitatively very different shapes of the interatomic potential (like 6-12 and screened 
Coulomb [HI [9]). In this sense, our results will be universal, although the molecular 
formation probabilities, to some limited extend, depend on the appearance of V{r). 
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To be specific we consider the pair collisions of and ^^Rb atoms with mass ration 
mi/m2 = 40/87 and the trapping frequency 002 = '^tt x 200 kHz for Rb. Hereafter we 
use the units fi = h = uq = 1 with cjq = 27r x 10 MHz. For the case of a decoupled 
CM motion, ui = uj2 = 0.02 in these units, whereas for the coupled case ui varies 
in the limits 002 < u^i < 2.2u2- Our focus is the analysis of the population dynamics 
of molecular resonance states (for reasons of brevity we will simply refer to them as 
molecular bound states with respect to V{r)) for ultracold atom-atom collision in the 
waveguide. 

3. Two-body bound states in a harmonic waveguide 

First we analyze how the molecular spectrum changes under the action of the confining 
waveguide (see Fig.l). 

We allow C12 to vary for fixed Cq = 1.847 ^20j and focus on the regime of the 
appearance of a bound state of V{r). In case of the absence of the confining potential 
the first bound molecular state Eq < appears at Cu — 0.13 and becomes increasingly 
bound with decreasing C12. To explore the dependence of the binding energies of these 
states on C12 in the presence of the waveguide we solve the corresponding eigenvalue 
problem of the four-dimensional Hamiltonian H{p^,r) employing the spectral 
method elaborated in ref. [21]. It is based on the computation of the autocorrelation 
function < il:{pR,r,t = 0)\'ip{pR,r,t) > where the initial state ip{pR,r,t = 0) can be 
considered as a test function of the spectrum. The solution ip{pji,r,t) of the time- 
dependent Schrodinger equation is obtained via our wave-packet propagation method 
[8]. For the case ui = UJ2 = uj, i.e. for a decoupled CM motion, every unconfined bound 
state with energy £:o(C'i2) transforms by the action of the waveguide into the spectrum 
£nin2{Ci2) which represents at Cu > 0.15 

BninziCu) (^^1 + ^^2) + 2?2iCJi + 2n2Co'2 = 2cj(l + Hi + n2) (6) 

the spectrum of two independent identical two-dimensional oscillators with radial 
quantum numbers Ui = 0, 1, ... and n2 = 0, 1, ... characterizing the transverse excitations 
of two noninteracting atoms in the waveguide. We use ni,n2 as labels of the molecular 
bound states eni,n2(C'i2) at fixed Cu correlating to the corresponding dissociation limit 
([6]) where they become exact quantum numbers. It is evident from ([6]) that for uJi = 002 
the excited states (ni,n2) become degenerate enin2 = ^n^nx- Since the coupling between 
the CM and the relative motion is absent the quantum numbers n and of the harmonic 
spectrum of the relative and the CM transverse eigenstates could be equally used to 
classify the spectrum. 

For L0\ 7^ 102 the presence of the coupling between the CM and the relative motion 
leads to a lifting of the above-mentioned degeneracy. In FiglHwe show the corresponding 
splitting £01 7^ ^10 of the first excited state and the shift of the ground state. 
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Figure 1. (color online) Molecular energy levels £0(^*12) and £ni,n2 (^"12) as a function 
of the parameter C12 of the atomic interaction in free space (open crossed circles) 
and in the waveguide. LOxjlfl = L02 = 0.02 for different transversal frequencies and 
LOi = UJ2 = 0.02 for the same frequencies. 



4. Molecular resonance states in waveguides 

The above-discussed spectral structure of the diatomic molecule in the waveguide 
permits us to analyze the dynamics of atomic pair collisions with respect to the formation 
of molecular final states. We choose herefore C12 = 0.109 and the energy of the colliding 
atoms to be e = ui + 002 + e\\ between the thresholds of the lowest ui + 002 and the 
first excited lo^ + 3ti;2 transverse channels (0 < < 2ti;2). Fig.l indicates that for the 
chosen potential and collision energies one can observe one bound state eoi of the closed 
channel ni = 0, n2 = 1 which becomes degenerate eoi = £^10 for cui = 002- Fig. 2 illustrates 
the time-evolution of the wave-packet in the course of the atomic collision for the case 
of CM nonseparability ojx 7^ 002 (Fig.2a) and CM separation ujx = uj2 (Fig. 2b). 

For uji 7^ U2 we observe that a considerable part of the scattered wave-packet 
(see the probability density distribution W{pR,r,t) at t = 8to) is located near r = 
after the collision and corresponds to a molecular bound state £01 which in the (n, A^)- 
representation is a mixture of the two states n = 0,N = 1 and n = 1,A^ = with 
dominating n = 0,N = 1 contribution. Note that during collision t ~ 3to the main 
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Figure 2. Time evolution of the probability density distribution averaged over the 
angular variables W{pii,r,t) — f \^{pii,r,9,(j),t)\'^{r'^ pii)~^ sm9d9d(j)- (a) For the 
waveguide with uji/2.2 = uj2 = 0.02 (b) for uji = u!2 = 0.02. For t = 8to corresponding 
to a time after the collision the insets show a more detailed view of W on the scales 
< W < 8x 10~^ corresponding to the continuum part of the wave-packet. e|| = 0.004, 
time is given in units of to — tt/uj2- 



part of the wave-packet is temporarily in the ground state eoo{n = N = 0), decays 
thereafter rapidly into the continuum but part of it goes into the excited molecular 
state £01- In contrast to this the case Ui = uj2 in Fig.2b (CM separation) shows an 
almost complete decay into the continuum: The remaining minor part near r = is 
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much smaller compared to the case ui ^002- 

To quantitatively investigate the molecular formation process accompanied by the 
CM excitation we calculate the population probability PAr(t) of the molecular bound 
states for = and = 1 

PN{t) = fj^ dr I dn \ dpRij{pR,r,t)^M{pR)? (7) 

where ^at are the two-dimensional oscillator states of the potential (l/2)(miCo'^ + 
m2CL'2)p|. and dVt = sin 6d6d(j). For the chosen interval of collision energies < ey < 2uj2 
and C12 = 0.109 the molecular states = and = 1 are the ones which dominate 
during the collision process (see Fig.l). Note, that the integration in ([7]) over the 
interatomic distance r is limited to the region r < = 10 of the action of the 6- 
12 potential. For sufficiently long times, i.e. after the collision process is concluded, 
all the unbound part of the scattered wave-packet resides outside this region. The 
complementary part is given by P/v(^) and can be interpreted as a molecular bound 
state. For — > cxd we obtain Po{t) + Pi{t) ~ 1 for any time and any energy within 
the above-provided interval. Fig.3 shows (cui 7^ UJ2) that during the time interval of 
close collision the wave-packet temporarily occupies the molecular ground state A^ = 
i.e. Pq becomes large. However, with further increasing time Pq decays to zero and Pi 
increases rapidly which corresponds to the population of the excited state A^ = 1. Pi 
decays very slowly for long enough times after the collision. We interprete this state 
as the first excited state ni = 0,n2 = 1 of the molecular spectrum (see Fig.l). For 
uji = UJ2 the coupling between the rzi = n2 = (A^ = 0) and ni = 0,^2 = 1(A^ = 1) 
molecular states is absent and the temporary population in the course of the collision 
of the molecular ground state decays back into the continuum after the collision (see 
behaviour of Po(^i = "^2)- In particular P\{u]x = UJ2) is negligible always. 

In FigJH we illustrate the mechanism of molecule formation for the ground state 
ni = n2 = 0{N = 0) with transversal Ui + UJ2 and relative longitudinal ey = /iWo/2 
energies. Under the action of the attractive tail —Cq/t^ of interatomic interaction the 
atoms accelerate up to the energy fiv^{t)/2 = e\\ + cui + 0^2 — £00 at the instant of 
collision. The occurence of the molecular state with energy ego (bound molecular state 
with respect to the ground transversal channel) as an intermediate in the course of the 
collision process can be seen by inspecting the intermediate probabihty density shown 
in Fig. 2a and Fig.3. The released energy subsequently transfers to the closed channel 
Ui = 0,n2 = 1{N = 1) via formation of a molecule in a CM excited state since this 
closed channel is coupled to the A^ = by the term W{pR, r) = piujI — uDrpR sin 6 cos 
(jl]) of the total Hamiltonian ([T]). This mechanism also explains the delay of the onset of 
the coUisional interaction observed in Figj3] in the absence of the coupling term uji = 002 
compared to the case ui 002- For ui = 002 the coUisional interaction happens later 
than for Ui 7^ U2 because the binding energy in the channel £:oo('^i 7^ ^2) is considerably 
lower than the binding energy in the channel £00 ("^i = 1^2) leading to a stronger atomic 
attraction and acceleration of the colliding atoms for ui ^ 002- A simple semiclassical 
estimate of the time difference between the onsets of coUisional interactions for the two 
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Figure 3. (color online) P/v(i) of the states A'^ = 0, 1 as a function of the collision 
time t. Full circles and squares correspond to ti;i/2.2 = uj2 — 0.02. Open circles and 
squares belong to oji = 0J2 = 0.02. ey = 0, 004, time is given in units of to = 'njioi- 



cases is in good agreement with the numerical result of Figl3J 

From the above scheme it also follows that if the energy release + tui + ti;2 — ^oo 
is equal to the excitation energy £01 ~ ^00 we have to expect a resonant enhancement of 
the molecular formation process i.e. we encounter the resonance condition 

e\\ + cui + u;2 = £01 (8) 

To analyze the energy dependence of the atom-molecule reaction we inspect the 
population probabilities P^[e\^ for the molecular states = and 1 as a function 
of the longitudinal energies e\\ (see FigjS]). The values P^is])^ = PN{£\\,t — > 00) were 
calculated after the collision at t = 8to, in the asymptotic region where Pi{e\\,t) is very 
slowly decaying if ui 7^ 002- Pi{£\\) exhibits a pronounced resonant behaviour whereas 
Po{£\\) shows a very weak energy dependence which is in agreement with the above- 
discussed reaction mechanism. The calculated position e\\ = Sr ^ 0.01 of the maximum 
of the population Pi{e\\) is shifted considerably compared to ErioJi 7^ uj-i) = £01(^1 7^ 
U2) — — ijj2 = 0.026 following from the resonance condition ([HD. We attribute this 
difference to the fact that we are using for the calculation of Pi{e\\) in eq.([7]) the oscillator 
wave function ^n{pr) instead of the exact transversal CM part of the molecular wave 
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Figure 4. Schematic illustration of the resonant formation of a diatomic molecule. 
Transfer of the released energy £|| + (uji +1^2)^ £00 to the excitation Eqi ^ £00 of the CM 
in the course of molecular formation. Solid lines indicate the positions of the molecular 
energy levels e„i^„2 fo^' the case uji/2.2 — uj2 — 0.02. The shifts of the levels due to the 
elimination of the coupling term at cji = cj2 = 0.02 are shown by dashed lines. 

function ni = 0,n2 = 1(A^ = 1). The resulting deviation from the exact resonance 
conditions ([8]) is particularly due to the missing effects of the CM coupling which 
provides the shift £01(^1 = ^2) — ^oii^^i 7^ ^^2) = 0.014 resulting in a significantly 
better agreement of the calculated position of the resonance with the resonant energy 
from eq.([H]). 

The molecular formation probability Pi depends also on Au = uji — uj2 and reaches 
a maximum at a value corresponding to the resonance condition ([8]). This probability 
can therefore be controlled by changing the transversal confinement frequencies for the 
different atomic species. By nature the process we observed is a molecular resonance in 
the atomic scattering and back decay of the molecule to the atom-atom continuum 
will occur. However, this resonance possesses a long lifetime due to its two-step 
character which is clearly demonstrated in FigOl Finally we note that Pi amounting to 
several percent should not obscure the fact that substantial molecular formation rates 
A = uaVqPi could be achieved by integrating over many collisions even for low linear 
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Figure 5. (color online) Pjv(e||) for the states iV = and 1 as a function of the 
longitudinal energy e|| before the collision. Full circles and squares correspond to 
cji/2.2 = UJ2 = 0.02. Open circles and squares correspond to uji = uj2 = 0.02. 



atomic densities n^. Thus, for ua = lO^cm ^ and T = lOOnK we obtain the estimate 
A ~ lOV^ 

5. Conclusions 

Perspectives concerning apphcations of the molecular formation mechanism are based on 
the fact that it represents a key ingredient for the preparation and following investigation 
of the dynamics of mixed atomic and dipolar molecular quantum gases in waveguides. 
The latter are expected to possess a novel collective excitation dynamics particularly 
in the ultracold regime discussed here for which only a few transversal channels are 
accessible. If wanted, one could stabilize the molecular gas and prevent back decay into 
atoms by a corresponding switch of the trapping fields such that the CM decouples from 
the interatomic motion finally. Open questions include the impact of inelastic atom- 
molecule and molecule-molecule collisions, which however goes beyond the scope of the 
present investigation. 
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